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Abstract 



The coherent states for a set of quadratic Hamiltonians in the trap regime are con- 
structed. A matrix technique which allows to identify directly the creation and annihilation 
operators will be presented. Then, the coherent states as simultaneous eigenstates of the 
annihilation operators will be derived, and they are going to be compared with those at- 
tained through the displacement operator method. The corresponding wave function will 
be found, and a general procedure for obtaining several mean values involving the canon- 
ical operators in these states will be described. The results will be illustrated through the 
asymmetric Penning trap. 
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m : 1 Introduction 

> 

^ ; As it was shown by Schrodinger in 1926 for the harmonic oscillator, the quasi-classical states 

are important for the description of physical systems in the classical limit (see e.g. flU). The 
catchy term coherent states (CS) was used for the first time by Glauber long after, when studying 
electromagnetic correlation functions J2l[3]|. With this application it was realized that the CS are 
useful as well in the intrinsically quantum domain. Indeed, the CS approach is nowadays widely 
employed for dealing with quantum physical systems. According to Glauber there are three 
equivalent ways to construct the CS for the harmonic oscillator. The first one is to define them 
as eigenstates of the annihilation operator. The second one is to build the CS through the action 
of a displacement operator onto the ground state. The third way is to consider them as quantum 
states having a minimum Heisenberg uncertainty relationship. These three properties can be 
used as definitions to build the CS for systems different from the harmonic oscillator. However, 
it is noteworthy that each of them leads to sets of CS which do not coincide in general ll4]-[T3l. 
In fact, even for the harmonic oscillator the third definition does not produce just the standard 
CS, since it includes as well the so-called squeezed states lfl4l[T5ll . 

In spite of its long term life, from time to time there are some advances which maintain 
the subject alive. This is the case, e.g., of the recently discovered coherent states for a charged 
particle inside an ideal Penning trap [16J. Since the corresponding Hamiltonian is quadratic in 
the position and momentum operators, one would expect that the CS appear as a generalized 
displacement operator acting onto the corresponding ground state. However, it is worth to notice 
that the Penning trap Hamiltonian does not have any ground state at all, since it is not a positively 
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defined operator. Despite, it was possible to implement in a simple way the corresponding CS 
construction. Thus, we need to take into account this Hamiltonian property when studying the 
coherent states of general systems. 

In this article we are going to address the CS construction for systems characterized by 
a certain set of quadratic Hamiltonians. The CS will be built up as simultaneous eigenstates 
of the corresponding annihilation operators, and also by applying a generalized displacement 
operator onto an appropriate extremal state. We will see that in the case of a positively defined 
Hamiltonian this extremal state will coincide with the ground state. In order to perform the CS 
construction, we need to find first the annihilation and creation operators. This task will be done 
by using a matrix technique, which generalizes the one employed in [fT6l (see also [fT7] - |20l ). In 
this way, we will simply and systematically identify the characteristic algebra of the involved 
Hamiltonians. Our procedure represents a generalization to dimensions greater than one of the 
standard technique to deal with the harmonic oscillator, which is closely related to the well 
known factorization method (see, e.g., H21 [12210 . 

The paper is organized as follows. In section 2 we will introduce a detailed recipe for 
systematically obtaining the annihilation and creation operators for quadratic Hamiltonians in 
the trap regime. The coherent states derivation shall be elaborated in section 3, while in section 
4 we will address the completeness of this set of CS, we shall obtain the mean values of several 
important physical quantities and the time evolution of these states. We are going to apply our 
general results to an asymmetric Penning trap in section 5, and our conclusions will be presented 
at section 6. 



2 Ladders operators for quadratic Hamiltonians 

Along this work we are going to consider a general set of n-dimensional quadratic Hamiltonians 
of the form 

H = ^rj T B V , (1) 

— » — » T — * -* 

where B is a 2n x 2n real constant symmetric matrix, r\ = ( X, P ) , and X, P are the 
n-dimensional coordinate and momentum operators in the Schrodinger picture satisfying the 
canonical commutation relationships [JQ, Pj] = iSij (notice that a system of units such that 
h = 1 will be used throughout this paper). The time evolution of the operator vector r](t) = 
W(t)r)U(t) in the Heisenberg picture is governed by 

^P- = tf(t)[iH, V ]U(t) = U\t)A V U(t) = A V (t), (2) 
at 

U(t) being the evolution operator of the system such that U (0) = 1, and A = JB, where J is 
the well known 2n x 2n matrix 

Mo). 

satisfying 

J T = -J, J 2 = -l 2n , det(J) = l, (4) 

in which l m represents the m x m identity matrix. The solution of Eq.© is given by 

v (t) = e At 7](0) = e A V (5) 
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In order to identify the annihilation and creation operators of if, we need to find the right 
and left eigenvectors of A. Since in general A is non-hermitian, its right and left eigenvectors 
are not necessarily adjoint to each other. 

Let us consider in the first place the 2n-th order characteristic polynomial of A, P(X) = 
det(A — A) = det(JB — A). Using Eqs.© and the fact that B is a symmetric matrix, we obtain 



P(A) = det(JB - A) = det[(JB - A) T ] = det(JB + A) = F(-A). 



(6) 



This means that if A is an eigenvalue of A, then — A also will be. Throughout this work we are 
going to denote the eigenvalues of A as A^ and — A&, taking A^ in the way 



Re(A fc ) > or Im(A fc ) > if Re(A fc ) = 0, k = 1, 



n. 



(7) 



Let us label as uj: and the right and left eigenvectors associated to the eigenvalues ±A* 
respectively, i.e., 



±X k u 



ff A = ±X h f, 



k ■ 



(8) 



Notice that both and fjj 1 can be determined from Eq.© up to arbitrary factors. Part of this 
arbitrariness will be eliminated by imposing two requirements. The first one is that the right 
and left eigenvectors be dual to each other, namely, 



(9) 



with j, k — 1, . . . , n and r, r' = +, — . The second condition, which is needed in order to recover 
the standard annihilation and creation operators for the one-dimensional harmonic oscillator, is 
to ask that the left eigenvectors involved in the commutators 



Ifk v, fkV] = ik, ik e C, 



are such that 



7* 



1. 



(10) 



(11) 



In this paper we are going to discuss just the case in which there are no degeneracies in the 
eigenvalues ±A fc so that the identity matrix l 2n can be expanded as (see Il23l ) 



n 



1-2 n 



(12) 



k=l 



with eg) representing tensor product. Then we get 



v(t) = 

n 

=£(■ 

k=l 

where L± = f±rj 



J2^®fk+ u k^fk 



k=l 



Xkt utLt + e~ Xkt ulL 



r]= y < 

k=i 



fkV + e 



-At* 



fkV] 



k^k )i 



(13) 
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It is worth to point out that, in the classical case, represent c-numbers which are related 
with the initial conditions. The time dependence of rj{t) is determined essentially by the 
values, which are complex in general. If Re(Afc) 7^ it can be seen that one of the two involved 
exponentials of the A;-th term in the previous relation diverges as t increases and, thus, the 
classical motion will be in general unbounded (see |[24"l - l27l0 . The only way in which this does 
not happen is that all the eigenvalues be purely imaginary so that the corresponding exponentials 
will just induce oscillations in time, and therefore in this case the classical evolution of the 
vector 77 (t) will remain always bounded. 

On the other hand, in the quantum regime the are linear operators in the canonical 
variables X, P. It is straightforward to show that their commutators with H reduce to 

[H,L±] = TihLt (14) 

In addition, it turns out that 

[LJ,Lj;] = [Lf,Lt] = Q, [LJ,L+] = 0, k ? 3. (15) 

However, 

[£fc>£fc]=7*^0, k = l,...,n. (16) 

Eq. (fT4j) implies that behave, at least formally, as ladders operators for the eigenvectors of H, 
changing its eigenvalues by =pAfc. However, this statement has to be managed carefully since 
it could happen that the action of onto an eigenvector of H produces something which does 
not belong to the domain of H, and in this case we would not get new eigenvectors of H. In the 
next section we will explore an interesting situation (the trap regime of our systems) for which 
the application of onto an eigenvector of H produce a new one with different eigenvalue. 

Let us point out that Eqs. (ll4H16l) imply that H can be expressed in a simple way in terms of 
{L± k = l,...,n} El. This is a consequence of the following general theorem: 

Theorem. If C is an irreducible algebra of operators generated by Lf which obey [Lf, Lf] = 
[LJ , L~] = 0, [L~, Lf] = ^idij with |7j| = l,i,j = 1, . . . , n, then an operator H G £ which 
fulfills relation (fT4l) can be written as 

where g G C. 

Demonstration. Actually, due to Eqs. (ll4H16l ) it turns out that 

commutes with for all k, thus with any function of them, and so g must be a c- number. □ 

From now on we are going to discard situations such that Re(Afc) 7^ for some k = 1, . . . ,n, 
restricting ourselves to cases in which the A& are purely imaginary for all k, i.e., we stick to the 
trap regime of our systems. 
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2.1 Algebraic structure of H in the trap regime 

Let us suppose that X k = iuk with oj k > 0, k = 1, . . . , n. Hence, \* k = — Afc, and since A is 
real, without loosing generality we can choose 

/* = (#)*, % = «)*■ (is) 

Since are linear combinations of the hermitian components of rj (Xi, . . . , X n , Pi, ... , P n ), 
it turns out that 

(LtV = L* (19) 

Note moreover that y k = y k , i.e., 7^ G M, and the use of Eq. (fTT|) implies that y k = ±1. 
Summarizing these results, Eq. (fl4)) becomes in this case 

[H,L±] = ±u k L±, (20) 



i.e., the modifications suffered by the eigenvalues of H through the action of the ladder operators 



are given by the real quantities ±uj k . In addition, H is factorized as (see Eq.(fT7l)) 



H = IkUkL+Ll + 00, 9o e R. (21) 
fe=i 

Notice that the previous summation either involves terms for which y k = 1, which are of the 
oscillator kind since they are positively defined, or terms with 7^ = — 1 which are of the anti- 
oscillator type since they are negatively defined. Thus, it is natural to define a global algebraic 
structure for our system (see, e.g., lfT6ll28~lQ . which is independent of the spectral details but 
has to do with the fact that any quadratic Hamiltonian in the trap regime can be expressed in 
terms of several independent oscillators, some of them indeed being anti-oscillators (compare 
Eq. (|2T]) N ). This global structure is characterized mathematically by identifying n sets of number, 
annihilation and creation operators of the system, {iV^, B k , B\], k — 1, . . . , n, in the way: 

Bl = L+, for 7fc = l, 

Bl = L~, for Tfc = -1, (22) 
k = 1, . . . , n, 

so that the standard commutation relationships are satisfied, 

[B v Bl\ =8 jk , [B J ,B k ] = [B},Bt] = 0, (23) 
[N k ,B k ] =-B k , [N k ,B\] = Bl j,k = l,...,n. 

Let us construct now a basis {\ni, . . . , n n ),rij = 0, 1, 2, ... , j = 1, 2, . . . , n} of common 
eigenstates of {N\, . . . , N n } (the Fock states) 

iV J -|ni,...,n T ,) = n J -|7ii,...,7i n ), j = l,...,n, (24) 

starting from an extremal state |0, . . . , 0), which is annihilated simultaneously by Bi, . . . , B n : 

B i |0,...,0) = 0, j = l,...,n. (25) 



B k 


= L~ k , 


B k 


= ^, 


N k 


= B\B k 
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If we assume that |0, . . . , 0) is normalized, it turns out that: 



ni,...,n n > = . (26) 



Moreover, Bj and £>j, j = 1, . . . , n, act onto \n\, . . . , n n ) in a standard way: 

Bj|ni, • • • ,nj-x,nj,n j+ x, . . . ,n n ) = y /n~\rii, . . . ,rij-x,Tij - l,n i+ i, . . . ,n n ), 

Bjlm, . . . ,nj_i,nj,n.,- + i, . . .,n n )= A /n J -fl|ni, . . . , n 3 -_i, njf 1, . . . ,n n ). (27) 

Now, in terms of the operators {Bj, Bj, j = 1, . . . , n} our Hamiltonian is expressed by 

n 

H = Y,lkUJ k BlB k + g' . (28) 



k=l 



It is clear that the Fock states \ni, . . . , n n ) are eigenstates of H with eigenvalues E ni ^ = 
7 1 t < j 1 r2 1 + . . . + ^y n u n n n + g' = ^(ni, . . . , ra n ). In particular, the extremal state |0, . . . , 0) 
has eigenvalue £7o,...,o = ffo- ^ n case mat 7fc = 1 for all k, then H — g' will be a positively 
defined operator, and the extremal state |0, . . . , 0) will become the ground state for our system, 
associated to the lowest eigenvalue £7o,...,o = g'o of H. On the other hand, if there is at least one 
index j for which jj = —1, then H — g' will not be positively defined, since the corresponding 
j-th term is of inverted oscillator type, and the state |0, . . . , 0) will not be a ground state for our 
system (however it keeps its extremal nature since it is always annihilated by the n operators 
Bj,j = 1, . . . ,n). 

Following [fT6ll28ll it is straightforward to see that, besides the global algebraic structure, 
there is an intrinsic algebraic structure for our system, characterized by the existing relationship 
between the Hamiltonian H and the n number operators N k : 

n 

H = E(N U ...,N n ) = J2 lkU k N k + g' . (29) 

k=l 

As in the examples discussed in lfT6ll28l . it turns out that this intrinsic algebraic structure is 
responsible for the specific spectrum of our Hamiltonian. On the other hand, the global alge- 
braic structure arises from the existence of the n independent oscillator modes for H, each one 
characterized by the standard generators {Nj, Bj, Bj}, j = 1, . . . , n. This global behavior al- 
lows us to identify in a natural way the extremal state |0, . . . , 0) = |0), which plays the role of 
a ground state although it does not necessarily has a minimum energy eigenvalue. Moreover, 
the very existence of the extremal state |0) is guaranteed by a theorem lT2"4l ensuring that if the 
operators {Bi, . . . , B n } obey the commutation relations given by Eq. (|23l) . then the system of 
partial differential equations 

(x\Bj\Q, ...,0) = (x\Bj\0) = 0, j = l,...,n, (30) 

has the square integrable solution 

(j) {x) = (f|0) = ce-fa**** = ce~^ TaS) , (31) 
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with c being a normalization factor. In this wave function, a = (a^ ) represents a symmetric 
matrix whose complex entries are found by solving the system of equations (1301) . leading to 

aaj = 4 j = l,...,n, (32) 

— # i 

where a,- and are obtained by expressing Bj and £>J as 

Bj =iP-dj + X ■ fy, B] = -id] -P + fit-X, j = 1, . . . , n. (33) 
The wave functions for the other Fock states can be found from Eq.d26l). 



3 Coherent States 

Once our Hamiltonian has been expressed appropriately in terms of annihilation and creation 
operators, we can develop a similar treatment as for the harmonic oscillator to build up the 
corresponding coherent states. Here we are going to construct them either as simultaneous 
eigenstates of the annihilation operators of the system or as the ones resulting from acting the 
global displacement operator onto the extremal state. 



3.1 Annihilation Operator Coherent States (AOCS) 

In the first place let us look for the annihilation operator coherent states (AOCS) as common 
eigenstates of the B/s: 

Bj\zi, . . . ,z n ) = z j \z 1 , . . . ,z n ), Zj E C, j = l,...,n. (34) 

Following a standard procedure, let us expand them in the basis {|ni, . . . , n n )}\ 

oo 

■ ■ ■ j Zn) ^ C-n\,...,n n ■ ■ ■ j Tin) ■ (35) 

ni,...,n Tl =0 



By imposing now that Eq. (|34l) is satisfied, the following recurrence relationships are obtained, 

Zj 

c ni,...,Ti,-,...,n„ = ;= c ni,...,n,'— l,...,n n j J ' = lj • • • i n i (36) 

which, when iterated, lead to 

n i 
Zj 

Cni,...,rij,...,n n = ...,0,...,n n i j 1, . . . , 71. (37) 



•nj\ 



Hence, it turns out that 



z^ 1 ... zy 



Cni,...,n„ — } , ==fC 0) ..o, (38) 

y/ml •■■n n \ 

where Co,...^ is to be found from the normalization condition. Thus the normalized AOCS 
become finally: 



/ 1 n \ oo 

\zx, . . . , z n ) = exp I -- 2j \zj\ 2 J Yl 



Z \ 1 • • • Z n n \ n li ■ ■ ■ J n n) 

Vnil ■■■n n \ 



up to a global phase factor. 
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3.2 Displacement Operator Coherent States (DOCS) 

The displacement operator for the j-th oscillator mode of the Hamiltonian reads 

D j {z j ) = exp(z j B}-? j B j y (40) 
By using the BCH formula it turns out that 

Dj(zj) = exp (~^p) ex P (zj B t) exp {-z*B 3 ) . (41) 

Now, the global displacement operator is given by: 

D{z) = D{ Zl , ...,z n )= D 1 {z 1 ) ■ ■ ■ D n {z n ), (42) 

where z denotes the complex variables z%, . . . , z n associated to the n oscillator modes. 

Let us obtain now the displacement operator coherent states (DOCS) |z) from applying 
D{z) onto the extremal state |0, . . . , 0) = |0): 

\ 3=1 / m,...,n n =0 v 



Notice that the AOCS and the DOCS are the same (compare Eqs.d39j) and (1431)). 



3.3 Coherent state wave functions 

In order to find the wave functions of the coherent states previously derived, we employ that 

3}-z*. 



[ Zj B] - z*Bj, z k B\ - z* k B k ] = V j, k. Thus: 



D(z) = exp(^ J B 1 t - zlBx) ■ ■ ■ exp(z n Bl - z*B n ) 

= exp[( Zl Bl + ■■■ + z n Bi) - (zjfli + • • • + z* n B n )} . (44) 



Using now Eq.([33l) we can write 

D(z) = e-*^-^) = e -|f.s e ffi.x e -if.p = e if.S e -ff.P e iS.X ) (45) 

where we have employed once again the BCH formula and we have taken 

f = 2Re[*J a 1 + ... + z*a n ], S = -2^*^ + . . . + z*J n \. (46) 
Now, it is straightforward to find the wave function for the coherent state |z), 

<f) z (x) = (x\z) = (x\D(z)\0) = e-* f ^^ s {x\e- ip - f \0). (47) 
Since the operator P is the coordinate displacement generator [Q]|, it turns out that 

(x\ e - ip - f = (x-f\, (48) 
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so that 

z (f) = e-^e^- s (x - f|0) = e-5 fg e ig -> (x - f). (49) 



A further calculation, using Eq.OTl). leads finally to 

z (x) = e -^ fTa+ig > f e( fTa+ig )->o(f). (50) 



Once again, it becomes evident that the extremal state is important in our treatment, since its 
wave function determines the corresponding wave function for any other CS. Moreover, as it can 
be seen from Eq.(|49l), the position probability density for the CS \z) becomes just a displaced 
version of the corresponding one for the extremal state |0). 

4 Mathematical and physical properties 

Let us derive next the completeness relationship for the previously derived coherent states. No- 
tice that, from the point of view of the analysis of states in the Hilbert space of the system, this 
is the most important property which our CS would have [|9l|TTl|29l[30l. This is the reason why 
several authors use it as the fourth coherent state definition, considering it as the fundamental 
one which will survive in time (see e.g. IfTTTl ). We are going to calculate as well some important 
physical quantities in these states. 

4.1 Completeness relationship 

A straightforward calculation leads to: 



|z)(z|d Z\ . . . d z n 



f \mu ; ..,m n )(nu...^ fr( 1 f z? ^ e -\ Zj \^ \ = ^ (51) 
n vmi!ni!---m n !n n ! - LJ - \7r J J / 

mi,ni,...,m n ,n n =0 n n ?=1 \ j / 



with 1 being the identity operator. Thus, the coherent states {|z)} form a complete set in the 
state space of the system (indeed they constitute an overcomplete set 0111321 ). This implies that 
any state can be expressed in terms of our coherent states, in particular, an arbitrary coherent 
state, 

l z ')=(^) J---J\z)(z\z')d 2 zi...d 2 z n , (52) 
where the reproducing kernel (z|z 7 ) is given by 



(z|z') = exp 



n 



2 

i=i 



(53) 



This means that, in general, our coherent states are not orthogonal to each other. Indeed, notice 
that inside our infinite set of coherent states only the extremal state of the system, |0) = |z = 
0) = |0, . . . , 0), is also an eigenstate of the Hamiltonian. 
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4.2 Mean values of some physical quantities in a CS 

Now we can calculate easily the mean values (X,) z = ( z \Xj\z), (fy)* = ( z \Pj\ z )>j = 1, • • • 
in a given coherent state |z) = \z±, . . . , z n ), as well as its mean square deviation in terms of the 
corresponding results for the extremal state |0). To do that, let us analyze first how the operators 
Xj, X? Pj, Pj are transformed under D(z). By using Eqs. (145 1461) it is straightforward to show 
that: 

£>t(z)X;D(z) = {X 3 + r,-) n , D\z)P?D(z) = {Pj + £,-) n , n = 1, 2, . . . (54) 
where we have used that, for an operator A which commutes with [A, B], it turns out that 

e A Be- A = B + [A, B] => e A B n e A = (B + [A, B]) n , n = l,2,.... (55) 

Thus, a straightforward calculation leads to 

(X j ) z = (z\Xj\z) = (Opt( z )X 7 -D(z)|0) = (X,-)o + IV (56) 

On the other hand, 

(X|) z = (X|)o + 2r,(X,) + r/. (57) 

Hence, 

(AX»2 = - (X,)^ = (X|) - (A^S = (AX»S. (58) 

Working in a similar way for Pj, it is obtained 

(P 3 ) z = (P 3 ) + £,-, (P, 2 ) z = (if ) + 2S i (P i ) + E/. (59) 

Then we have as well that 

(AP,) 2 = (AP^, (60) 

i.e., the mean square deviations of Xj and Pj in the CS |z) are independent of z. 

In order to end up this calculation, the mean values (Xj) Q , (P/)o, (Xpo, and (P/)o for 
j = 1, . . . , n are required. Let us describe now the procedure to find these An quantities. The 
first 2n, (Xj) , (Pj)o, can be easily found by recalling the definitions of B k , B\ (see section 2. 1) 
and using the fact that their mean values in the extremal state 1 0) always vanish for k = 1, . . . , n: 

(S fe )o = <P4)o = 0, fc = l,...,n. (61) 
This is equivalent to the following linear system of 2n homogeneous equations 

fk(v)o = ti(v)o = 0, k = l,...,n. (62) 

Since the left eigenvectors , k = 1, . . . , n, are linearly independent, the only solution for the 
2n unknowns (?/) is the trivial one, i.e., (Xj) = (Pj)o = 0. It is worth to point out that this 
result simplifies Eqs. (|56l57l59l) . 

On the other hand, the mean values of the quadratic operators Xf, Pf, i — 1, . . . , n, in the 
extremal state |0) can be obtained from evaluating the corresponding quantities for the several 
non-equivalent products of pairs of annihilation Bj and creation B\ operators. It is important to 
mention that these products should have the appropriate order to use the fact that Bj annihilates 
|0) and B\ annihilates (0| (if the product involves one Bj it should be placed to the right while 
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if it involves one B* k it should be placed to the left). In general, we will get n{2n + 1) non- 
equivalent products of pairs of operators Bi, Bj, i,j = 1, . . . ,n: n{n + l)/2 products of kind 
BiBj,j = l,...,n, i < j; n(n + l)/2 products of kind BjBLj = l,...,n, % < j; n 2 
products of kind B}Bj, i,j = 1, . . . , n. The mean values in the extremal state |0) will lead to 
an inhomogeneous systems of n{2n + 1) equations, with the same number of unknowns. When 
solving this system we will get (X 2 ) , {Pf)o, j = 1, • • • , n, and the mean value of any other 
product of two canonical operators Xi, Pj. 

It is customary nowadays to group the mean values of the quadratic products of the operators 
Xi, Pj in the coherent state |z) in a 2n x 2n real symmetric matrix a(z), called covariance 



matrix, whose elements are given by (remember that r] = (X%, . . . , X n , Pi, . . . , P n ) T )'- 

aij(z) = + %r?i) z - (r]i)z(r]j)z, i,j = l,... 2n. (63) 

A straightforward calculation leads to 

o-y(z) = cry(O) = + r)jr)i) = crij, (64) 



where we have used that (rji) o = 0, i — 1, . . . , 2n. The conclusion is that the covariance matrix 
in our coherent state |z) is once again independent from z and depends just of the extremal state 
|0). Notice that the number of independent matrix elements cr^ (n(2n + 1)) coincides with the 
number of unknowns which are determined from the set of n(2n + 1) independent equations 
associated to the mean values of the quadratic products of B k in the extremal state |0). 

Once the covariance matrix is determined, the generalized uncertainty relation can be eval- 
uated J3M1 

@ii ^n+in+i &in+i — 45 * \, . . . ,Tl, (65) 

which coincides with the Robertson-Schrodinger uncertainty relation (see e.g. 11331 ). 

Let us end up this section by calculating the mean value of the Hamiltonian in a given 
coherent state |z). Equation (|28l ) leads to 

n 

(H) z = (z\H\z) = ^'j k uj k \z k \ 2 + g' . (66) 

k=l 

In order to get (H 2 ) z , let us notice that 

n n n 

h2 =Y1 l 3 lkOJ 3 uJ k B)B\BjB k + u 2 k B ] k B k + 2g' ^ lk u k B ] k B k + g' \ (67) 

j,k=l k=l k=l 

Thus we get 

n n n 

(H 2 ) z = Y ljlk^j^k\zj\ 2 \z k \ 2 + ^uj 2 k \z k \ 2 + 2g' ^2'j k uj k \z k \ 2 + g' Q 2 . (68) 

j,k=l k=l k=l 

Hence, 

n 

(AH) 2 z = J2"l\z k \ 2 - (69) 

k=l 

Notice that, for one-dimensional systems (n = 1), this expression reduces to the standard 
one for the harmonic oscillator (see e.g. [HI). 
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4.3 Time evolution of the CS 



Suppose that at t = our system is in a coherent state |z) . Thus, at a later time t > the evolved 
state is found by acting on |z) with the evolution operator of the system U(t) = exp(-iHt). 
By making use of Eq.(|29l) it turns out that 

n 

U(t) = e -^o*TJ e - i7 ^ fcArfe . 

k=l 

Hence, 

U(t)\z)=e~^ t \z 1 (t),...,z n (t)) = e-< t \z(t)), (70) 

where Zj(t) = e~ l/yjU>jt Zj = \zj\e % ^ 3 ~ l3U)jt \ Equation (|7Q|) implies that a coherent state |z) 
evolves in time into a new coherent state |z(i)) = \zi(t), • • • , z n (t)), where the j-th degree of 
freedom Zj(t) just rotates at its characteristic frequency Uj (clockwise if jj = 1 and counter- 
clockwise if jj = —1). 

4.4 Gazeau-Klauder coherent states 

At this point, it would be interesting to check if our CS belong to the class introduced recently 
by Gazeau and Klauder [|36ll . Using their notation, for a system with a Hamiltonian % such that 
the ground state energy is zero, the Gazeau-Klauder CS {| J, 9), J > 0, —00 < 9 < 00} obey 
the following properties: 

(a) Continuity: (J', 9') -> (J, 9) \ J',9') -)■ |J,0). 

(b) Resolution of unity: 1 = f\J,9) (J, 9\d/i(J, 9). 

(c) Temporal stability: e~ lHt \ J, #) = \J,9 + tut), u = constant. 

(d) Action identity: ( J, J, 9) = uJ. 

Concerning the first property, it is straightforward to check that our coherent states given in 
Eq.(l43l) are such that |z') — > |z) as z' — > z, i.e., they are continuous in z. As for the second and 
third properties, both were explicitly proven in sections 4.1 and 4.3 respectively. It remains just 
to analyze if it is valid the action identity given in (d). Let us notice first of all that it is valid for 
each partial Hamiltonian H k = r y k ujkN k of our system, 

(z\H k \z) = 7fcWfckfc| 2 , 

which is time-independent. Therefore, property (d) becomes valid for each degree of freedom 
separately and thus it is valid for our global system with the natural identification J k = \z k \ 2 , 
9 k = arg(^ fc ) so that 

n 
k=l 

We conclude that our CS of Eq.(l43l) become as well an n-dimensional generalization of the 
Gazeau-Klauder CS if we express each complex component z k of z in its polar form (the polar 
coordinates essentially coincide with the canonical action-angle variables for the corresponding 
classical system). 
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5 Asymmetric Penning trap coherent states 



Let us apply now the previous technique to the asymmetric Penning trap. Such an arrangement 
can be used to control some quantum mechanical phenomena 11371 as well as to perform high- 
precision measurements of fundamental properties of particles. Moreover, it is a quite natural 
system to analyze the decoherence taking place due to the unavoidable interaction of the system 
with its environment [|38ll3~9ll . Since the asymmetric Penning trap becomes the ideal one when 
the asymmetry parameter vanishes ||40l - l42ll . it will be straightforward to compare these results 
with those recently obtained for the ideal Penning trap |fT6l (see also Il33ll4"3"1i 44|). 

The Hamiltonian of a charged particle with mass m and charge q in an asymmetric Penning 
trap reads 

H = £ + Y (XPy ~ YPx) + T (w ** 2 + ^ + ulz% (71) 

u c = qB/m and u z being the cyclotron and axial frequencies respectively, and the frequencies 
oj x , u y are given by 



4 



2 



1 



(72) 



where |e| < 1 is the real asymmetry parameter and we are denoting P = (P x , P y , P Z ) T , X = 
(X, Y, Z) T . Without loosing generality IfToll . from now on we will assume that m = 1. 

As it was seen at section 2, the main role in our treatment is played by the matrix A such 
that [iH, r]] = Arj. We choose here i] = (X, Y, P x , P y , Z, P Z ) T so that 



A 



The eigenvalues (A) of A are 

Ai 



f 


-w c /2 


1 








o\ 


u c /2 








1 








-< 








-Wc/2 











-cu 2 

V 


u c /2 


























1 


V o 













0/ 



100, 



A, 



-YV2-5 + R = ioj x 



A2 = ^V2 - 5 - R = iu 2 , 
R = y/A(l -5) + 5 2 e 2 , < 5 = < 1, 



(73) 



(74) 



and their corresponding complex conjugate. The right (u) and left (/) eigenvectors of A be- 
come 

i 2 + 5e + R iu c 2(1 - 5) - 5e + R , n ^ x 1 

U{ — S\ I — — , , ; , 1, 0, 



U 9 



u. 



S2 



uj c (5e + R) wi 5e + R 2u x 5e + R 

4 i 2 + 5e-R iu c 2(1 — 8) — 8e — R 

uo c (8e — R) -1 0J2 8e — R ' 2uj2 8e — R 

T 



,1,0,0 



s 3 0,0,0,0, ,1 
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/+ = h {^{R - 5e), 02(1 -S)+5e + R],-^-(2-6e + R), 1, 0, o) , 

ft = t 2 (j(-R - 6e), ^[2(1 -S) + Se- R], -|^(2 - fe - it!), 1, 0, o) , 

f+ = t 3 (0,0,0,0,^,1), (75) 

where Sy,tj G C,j = 1, 2, 3. The requirement that the right and left eigenvectors be dual to 
each other implies 

1 / 5e\ 1 / 5e\ 1 

On the other hand, up to some phase factors, the condition imposed by Eqs. dlOll II) leads to, 

1 1 1 

h ~ ^2i{f la f lc = J^ ) ' h ~ ^2i(f 2b -ftjS J ' 3 ~V^T (77) 

where we are denoting /{•" = h(f la , fib, fu, 1, 0, 0), / 2 + = t 2 (f 2a , fa, he, 1, 0, 0) in order to 
simplify the notation (compare Eq.(f75T)V Moreover, the crucial signs for us to conclude that our 
asymmetric Penning trap Hamiltonian is not positively defined become: 

7i = !> 72 = -1, 73 = 1- (78) 
Thus, our annihilation operators take the form (see Eq. (122)0 : 



B X = L\ = tx(f la X-f lb Y-f lc P x + P y ), 

B 2 = L+ = t 2 (f 2a X + f 2b Y + f 2c P x + P y ), (79) 
5 3 = L 3 = t 3 (-iu 3 Z + P z ). 

From these operators and their hermitian conjugates, it is straightforward to identify the otj and 
f3j, j = 1, 2, 3 which allow us to find the matrix a such that actj = j3j. Its matrix elements ay- 
become now: 

. fla — /2a . fib + f2b . /lc/26 — J 2c fib 

a n = - % - f — 7~c~-< a i2 = a 2i = i~ f — — j-, a 22 = I — , 

Jlc + J2c Jlc + he he + he 

a 33 = u 3 , ai 3 = a 3i = a 23 = a 32 = 0. (80) 

It can be shown that an, a 22 , a 33 G IR + while a i2 is purely imaginary. Thus the extremal state 
wave function of Eq.dSTT) acquires the form: 

O (^) = cexp ^-^a n x 2 - ^a 22 ?/ 2 - a 12 xyj exp ^-^a 33 z 2 ^ . (81) 

The associated eigenvalue becomes £7o,o,o — (^1 — ^2 + w 3 )/2. 

Concerning the coherent states z 2 > 23), the general treatment developed in section 3 is 
straightforwardly applicable, and their explicit expressions are given by Eq.(l43l) with n — 3. 
Their corresponding wave functions are given by 



(x 



(x\z) = e-^/V 2 -*^ - r ls y - r 2 , z - r 3 ), (82) 
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where 
f = 2 



iti/i c Re [zi] - it 2 f 2c Re [z 2 ] 
-iilm [zx] - t 2 Im [z 2 ] 
-t 3 lm [z 3 ] 



hf la Jm [zi] + t 2 f 2a lm [22] 
-iti/i&Re [^i] + it 2 f 2b Re [z 2 ] 
t 3 u 3 Re [z 3 ] 



(83) 



The mean values (Xj) z , (Pj) z , immediately follow from Eqs. (l56[|59l) with (Xj)n = (Pj)o 
0, i.e., 



31 



J 



1,2,3. 



(84) 



As for the mean values of the quadratic operators in the extremal state, we have solved the 
system of equations arising from the null mean values of the products of pairs of annihilation 
Bj and creation B>\ operators. We get 



and the crossed products 



(XP X ) - -, 



2a n ' 
1 

2a 22 
1 



(Px)o — o ( a H ~ 



a 22 

^33, 



a 22 



"12 
an 



(85) 



Wo = p, 
z a 22 

(ZP X ) = 0, 



(YPy) 



2 • 



(ZPy) = 0, 



(yP,)o 
(ZP g ) 



0. 

0. 

i 

2' 



(86) 



Therefore, using Eqs. (l58ll60l) we get the Heisenberg uncertainty relationships 



(AX)^(AP^ 
(^Z)l(AP z )l 



(^Y)l(AP y )l 
1 

4' 



1 + 



|Ql2| 

ana 22 



> 



(87) 



while Eq.(|69l) with n = 3 gives the mean square deviation for the Hamiltonian. 

Once we have calculated the mean values of the quadratic products given in Eqs. (l85l86l ), 
it is straightforward to evaluate the covariance matrix elements of Eq.(|64l). With the ordering 

77 = (X, Y, P x , P y , Z, P Z ) T , it is obtained: 



a 



/(AX) 2 








iai2 
2an 








\ 





(AF) 2 


io-12 
2a22 
















2a22 


(AP,)S 













2an 








(AP,) 2 






















(AZ) 2 





J 


V 














(AP Z 



(88) 



15 



Notice that this covariance matrix is non-diagonal. However, since a ri = ct 2 4 = °"56 = 0, it 
turns out that the generalized uncertainty relations of Eq.(l65l) reduce to the Heisenberg uncer- 
tainty relations given in Eq.(l87l). 

A plot of ( AX) Z ( AP X ) X as a function of the parameters e and 5 is given in Fig. 1 . As it can be 
seen from Eqs. (|8QI87l) and from Fig.l, the coherent states minimize the Heisenberg uncertainty 
relationship for e — 0, which coincides with the results recently obtained for the ideal Penning 
trap [fT6l . However, for e ^ it turns out that (AX) Z (AP X ) Z > 1/2. Notice that the same plot 
will appear for the uncertainty product (AY) z (AP y ) z . 




i.o - L0 



Figure 1: Heisenberg uncertainty relationship (AX) Z (AP X ) Z for the asymmetric Penning trap coherent states as 
function of the real dimensionless parameters |e| < 1, < S < 1. 



6 Concluding remarks 

In this work we have proposed a systematic technique to find the CS for systems governed by 
quadratic Hamiltonians in the trap regime. To do this, we introduced a prescription to identify 
in a simple way the appropriate ladder operators which play the same role as the annihilation 
and creation operators for the 1-dimensional harmonic oscillator. These operators allowed us 
to generate the eigenvectors and eigenvalues for the Hamiltonian departing from the extremal 
state, the analogue of the ground state although it is not necessarily an eigenstate associated to 
the lowest possible eigenvalue. The explicit expression for the extremal state wave function was 
as well explicitly calculated. 

For systems governed by this kind of Hamiltonians the two algebraic CS definitions (ei- 
ther as simultaneous eigenstates of the annihilation operators or as resulting from the action of 
the displacement operator onto the extremal state) lead to the same set of states. The explicit 
expression for the corresponding wave functions has been also derived. 

We have calculated explicitly the mean values of the position and momentum operators 
in an arbitrary coherent state. Moreover, we have provided as well a prescription to obtain 
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algebraically, by solving a linear systems of equations, the mean values of the quadratic products 
of these operators in the CS. 

Through this method we have found the asymmetric Penning trap coherent states and we 
have explored some of their physical properties. In particular, it is worth to point out that, in 
general, they do not minimize the Heisenberg uncertainty relationship. The differences from 
the minimum are induced by the deviations of the axial symmetry which the ideal Penning trap 
has (measured by the asymmetry parameter e). 
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